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Abstract. We study the extent to which spin assignments of new particles produced at the LHC can be 
deduced in the decay of a scalar or fermion C into a new stable (or quasi-stable) particle A through the 
chain C B^q,B^ AW^ £^ve where £ = e,ii. All possible spin assignments of the particles 
A and B are considered. Explicit invariant mass distributions of the quark and lepton are given for each 
set of spins, valid for all masses. We also construct the asymmetry between the chains with a W~ and 
those with a . The KuUback-Leibler distance between the distributions is then calculated to give a 
quantitative measure of our ability to distinguish the different spin assignments. 



1 Introduction 

While the Standard Model (SM) has been remarkably suc- 
cessful to date, new physics is expected around the TeV 
scale, for example to cancel the large contributions to 
the Higgs mass thereby solving the Hierarchy problem. 
Whatever form this new physics takes, we expect to find 
new particles. The issue of deducing the spin of these new 
particles from experimental data has become increasingly 
important with the rise in popularity of supersymmetric 
(SUSY) extensions to the SM. These models assign to SM 
partners different spin to that of the corresponding SM 
particle. 

This is in contrast to another possible SM extension, 
Universal Extra Dimensions (UED) T] where each SM 
partner has the same spin as its SM counterpart. In these 
models all fields propagate into at least one extra dimen- 
sion forming Kaluza-Klein towers of new particles with in- 
creasing mass but otherwise identical quantum numbers. 
From this construction, a typical UED mass spectrum is 
very degenerate. There are also other possible extensions 
to the Standard Model such as Little Higgs Models [2] 
where the Higgs field is a pseudo Nambu Goldstone boson 
from a broken symmetry group. These models often fea- 
ture different spin assignments to new particles, such as 
new scalars without a direct SM counterpart. 

Often, studies of spin are considered in the context of 
a linear electron collider. However, Barr [F^ (see also [Al) 
showed that it was possible to deduce such information 
at the Large Hadron Collider (LHC). He demonstrated 
that one could distinguish between the case where parti- 
cles had SUSY spin allocations and where the particles 
were all effectively spinless. This work was extended in [5l 
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l6l[7ll8l[9lfT0] to demonstrate that spin studies were a use- 
ful tool to distinguish between SUSY and UED. (It was 
first pointed out in [TT] that these two models could mimic 
each other.) Recently [T2l, the technique was extended to 
cover all possible spin assignments in the cascade decay 
of a quark partner via opposite-sign-same-flavour (OSSF) 
leptons. This had much wider applications as it was no 
longer constrained to spin effects only in the Minimal Su- 
persymmetric Standard Model (MSSM) and UED. A sim- 
ilar study had previously been applied in [13] to the pair 
production of top quark partners each decaying straight 
to a top and a stable particle. 

These studies have concentrated on the quark partner 
cascade decay (or gluino decay leading to this), top part- 
ner production and Drell-Yan production of lepton pairs 
and their subsequent decay. Here we study the electroweak 
decay of a quark partner via a W boson decaying lepton- 
ically. In the MSSM, this decay chain often has a higher 
branching ratio than the cascade decay via a X2 which is 
more frequently studied. In [T4l, it is suggested that this 
could be the most promising channel for spin discrimina- 
tion. Here, we consider all possible spin assignments so as 
not to constrain ourselves to a particular model. We as- 
sume that these chains have been identified, and that the 
masses of the particles involved are known. The spin corre- 
lations in the chain depend on the charge of the W boson, 
so we consider the two charge assignments separately. 

In section [21 we discuss all possible spin assignments 
in the decay chain and the resulting matrix elements. 
In section [3l spin correlations are discussed in terms of 
the invariant mass distributions of the quark and lepton. 
The full analytical formulae valid for any mass spectrum 
are calculated. We then form an asymmetry between the 
chains with a W~ and the chains with a W^. These distri- 
butions are plotted and discussed. In section [H we quan- 
tify the results of the previous section using the KuUback- 
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Leibler distance method introduced in [12 . This gives a 
lower limit on the number of events required to discrim- 
inate between any two of the spin allocations at a given 
level of confidence. These lower limits do not include back- 
ground or detector effects as these will vary between ex- 
periments and we wish to remain general. They do how- 
ever provide a handle on the feasibility of distinguishing 
two particular curves. This analysis is applied to the ob- 
servable processes individually, and then combined. The 
conclusions are in section [5] before the more lengthy for- 
mulae and a discussion of higher derivative vertices which 
are in the appendix. 



2 Spin assignments 

We will consider the decay of a heavy colour-triplet scalar 
or fermion C of the form C B^q, AW^, 
i^i^i (figure [T]), where £ = e,^. Chains like this can occur 
in SUSY, UED or the Littlest Higgs model with T-parity. 




Fig. 1. The decay chain under consideration. 

We will assume that particle A is a stable or long-lived 
heavy massive new particle and that the masses of the new 
heavy particles A, B and C have all been measured. All 
possible spin configurations are listed in Table [l] together 
with the labels which will be used in the rest of the paper. 



Label 


C 


B 


A 


SFF 


Scalar 


Fermion 


Fermion 


FSS 


Fermion 


Scalar 


Scalar 


FSV 


Fermion 


Scalar 


Vector 


FVS 


Fermion 


Vector 


Scalar 


FVV 


Fermion 


Vector 


Vector 



Table 1. Possible spin configurations in the decay chain (figure 



The SFF chain corresponds to SUSY spin assignments 
while FVV corresponds to the spin assignments in a UED 
model, or a Littlest Higgs model with T-parity [TS]. The 
other spin assignments correspond to non-minimal ver- 
sions of these or other models. The UED masses derived 
from [inj do not allow a decay chain of this form to pro- 
ceed for values of the compactification radius accessible 
at the LHC, however, these were calculated under the as- 
sumption that the orbifold boundary kinetic terms vanish 
at the cut-off scale. This is not necessarily the case and for 
different values of these parameters it is possible that the 



decay chain would still proceed. Indeed, it is the freedom 
in choice of parameters which makes it unlikely that these 
models could be distinguished by studying mass spectra 
alone. 

It is necessary to make some assumptions about the 
structure of the vertices in the chains, except for the SFF 
chain where these are well-determined in the MSSM. We 
are not concerned with overall numerical factors as the 
distributions are normalised to integrate to 1. When we 
consider the FSS, FSV, FVS and FVV chains, the B-W- 
A vertex structure is uniquely determined if we do not 
consider higher dimensional couplings like those induced 
from loops. In the FSS chain, it is of the form (p — q)^ 
where p and q are the incoming momenta of the scalars; 
in the FSV and FVS chains it is of the form g'^'^, where 
fj, is the index corresponding to the W and v is the index 
corresponding to the other vector particle {A in FSV or B 
in FVS), while in the FVV chain the triple vector vertex 
takes theiormg''''{pi-p2y+g^P{p2-P3T+9'"'{P3-PiT- 
These are the structures considered here — a discussion 
of possible alternatives is in appendix [Bl 

The structure of the C-q-B vertex is not so well deter- 
mined in these chains and in principle contains a factor of 
(1 + 075) where a is an arbitrary constant. In the massless 
q limit, the final distributions are in fact independent of 
a for the FSS and FSV chains, but this is not the case 
for FVS and FVV. For these chains, where necessary, the 
constant a value has been taken to be —1, thereby forc- 
ing the q to be left-handed. This value is justified because 
most models beyond the standard model have two excita- 
tions for each fermion — one coupling to the left-handed 
fermion and one coupling to the right. As they have the 
SM as a low energy limit, it is usually the one associated 
with the left-handed fermion which undergoes decays of 
the type in figure [1] (especially for the light quarks we 
consider where left-right mixing is expected to be small). 
In particular, the FVV chain has the spin structure found 
in UED where this is the case. 



3 Spin correlations 

In the chain, there are only two observable emitted par- 
ticles, the quark (jet) and the charged lepton. This gives 
one observable invariant mass-squared: m^^ = [pq -\- pi)^ . 
We define the angle 9 to be the angle between the quark 
and A in the rest frame of i?, and ■(/' to be the angle be- 
tween the lepton and A in the rest frame of the . We 
then define (p to be the angle between these two planes. 
Then, 

ml, = ^m|(l-X)((l + y-Z)(l-cos0cosV) 

+ ^{i + Y - zy - 4y(cos e ~ cos V') 

—2^/Y sin 61 sin -0 cos 0^ , (1) 

where the mass-squared ratios X^Y, Z are X — ra^/m^, 
Y = m^/m^g and Z = m\/m^g. These must satisfy ^/Y + 
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< 1 by energy conservation and so the quantity in the 
square root is always non-negative. The maximum value of 
^le is Jx^l{^-X){{l + Y-Z) + ^{l + Y ^ Z)2 -4Y) 
which occurs when {0,ip) = (0, tt). 

In order to keep a manageable expression, we define 
the scaled invariant mass as 

which lies in the interval 

0, 2((1 + Y-Z) + ^{1 + Y - Zy -4Y) . 

The analytical expressions valid for any particle masses 
are discussed in appendix]^ however, in order to plot the 
functions we must choose values for the masses of A, B 
and C . If we consider this chain in a SUSY scenario, we 
have the particle assignments given in table [2l The masses 
given are those at the Snowmass Benchmark points, SPS 
la, SPS 2 and SPS 9 [E]. SPS la and SPS 2 were chosen 
as the points with the biggest difference in their spectrum, 
while the AMSB point, SPS 9, was chosen as an example 
of a heavier chargino which allows for a greater difference 
between the mass ratios Y and Z . 





C 

UL 


B 

xt 


A 

x1 


SPS la 


537 


378 


96 


SPS 2 


1533 


269 


79 


SPS 9 


1237 


876 


175 



Table 2. The mass spectra (in GeV) considered in this paper. 



The spin correlations in the chain where the quark 
partner decays through a, W~ has different spin correla- 
tions to that in which the quark partner decays through 
a as one has a charged lepton, the other a charged 
anti-lepton (this sends a ^ —a). This means we have two 
processes to consider: 

Process 1: {q,W^} = or {m,W^+} 

Process 2: {q,W} = {d,W+} or {d,W-}. 

Here u stands for either an up or a charm quark and d 
stands for a down or a strange quark. We do not include 
bottom and top quarks since b and t final states should 
be distinguishable from those due to the lighter quarks. 
We may then work in the massless approximation. For 
the FSS and FSV chains, processes 1 and 2 give the same 
distribution as the scalar does not carry spin information 
down the chain. 

Figure [2] shows the invariant mass-squared distribu- 
tions for both processes for the SPS la mass spectrum, for 
the five spin assignments given in Table [TJ Here, the dis- 
tributions are plotted as dP/dm throughout, as opposed 
to dP/dm^ as was done in 12], as the phase space is not 
flat in any such simple mapping of the invariant mass. The 
phase space curve (the case where all particles are treated 
as spinless) also depends on the masses in the chain, and 
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Fig. 2. Invariant mass distributions for SPS la: Process 1 (top) 
and Process 2 (bottom). 



SO is indicated on the dPi/dm plot for each mass spec- 
trum (marked No Spin) . Figures [3] and |4] show the same 
distributions for the five spin assignments, for the SPS 2 
and SPS 9 mass spectra. 

These plots show that all the curves have a similar 
overall shape, but with some differences due to the dif- 
ferent Lorentz structure. The exact effect can be seen in 
the equations in appendix|3 However, quantitative state- 
ments can be made. For example, the SFF (MSSM) curve 
peaks slightly to the left (right) of the others in Process 
1 (2) for all mass spectra, although to different degrees. 
Also, the FSS and FSV curves are very similar particularly 
at SPS la and 9. 

From these curves for processes 1 and 2, we can con- 
struct the distribution of processes through a W~ and the 
distribution of processes through a , which are the dis- 
tributions which would actually be observed. If an l~ is 
observed, the chain must have started with the partner 
of either a down-type quark, or an up-type antiquark. We 
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Fig. 3. Invariant mass distributions for SPS 2: Process 1 (top) Fig. 4. Invariant mass distributions for SPS 9: Process 1 (top) 
and Process 2 (bottom). and Process 2 (bottom). 



define r^. = 1 — r^* to be the fraction of chains with an £^ 
which begin with the partner of a down-type quark. Sim- 
ilarly, we define r„. = 1 — rj. to be the fraction of chains 
with an £^ which begin with the partner of an up-type 
quark. The q£^ distributions, dPzp/dm, are given by: 



dp_ 


dPi 


dPa 


dm 


dm 


dm 


dP+ 


dP2 




dm 


" dm 


dm 



(3) 



No distinction between flavours of quarks was required 
in the earlier studies [51l51[71fT^ of the cascade decay of a 
quark partner. In these the quark partner decayed straight 
into a quark and a neutral particle so no charge informa- 
tion of the original quark was transmitted to the rest of 
the chain making the results flavour independent. 

The MSSM scenarios in table imply the values of 
the fractions in table [3] at the LHC, i.e. in pp collisions at 



14 TeVQ Therefore these values represent studying models 
with the MSSM flavour structure but different spin assign- 
ments. We see that at SPS 9, the effect of having more up 
quarks than down quarks in the proton is dwarfed by the 
latter's larger branching ratio to a chargino. This is caused 
by the large value of the MSSM parameter fi enhancing 
the effect of large tan/3 in the chargino mixing matrices. 
The resulting plots are shown in figures [5] - [71 

We see that at SPS 9 the plots are nearly identical 
due to the extreme values of r^* and rj* there. There is 
greater variation in the individual curves at SPS 9 than 
for the other two mass spectra. Our ability to distinguish 
the curves is discussed in section |4l 



^ These results were obtained from the Herwig [T8][T9][20] 
Monte Carlo at parton level, corresponding to the leading- 
order QCD production processes and MRST parton distribu- 
tion functions 121] . They are not sensitive to details of the 
Monte Carlo, higher-order corrections or PDF uncertainties. 
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Spectrum 




rw 






SPS la 


0.860 


0.140 


0.469 


0.531 


SPS 2 


0.900 


0.100 


0.911 


0.089 


SPS 9 


0.998 


0.002 


0.072 


0.928 



Table 3. Numerical calculation of fractions using Herwig. 




Fig. 5. Observable invariant mass distributions for SPS la: 
P- (top) and P+ (bottom), see equation (|3}. 



We combine the information from both chains together 
by forming the asymmetry of the normahsed distributions 
given by 

dP _ dP+ 

^ - dP_ dP+ ■ V4j 
dm^ dm-^ 

The resulting plots are given in figure [51 

The asymmetry at SPS la (figure [U top) shows a dif- 
ference in the behaviour of the SFF, FVS and FVV curves. 
With a 10% level of asymmetry, we can be optimistic 
about distinguishing the SFF (MSSM) curve. The differ- 
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Fig. 6. Observable invariant mass distributions for SPS 2: P_ 
(top) and P+ (bottom)), see equation Q. 

ences at very high and low m cannot usually be used as 
this is where experimental statistics are often much worse. 

In the asymmetry plot for SPS 2 (figure O middle), 
it is unlikely that we would be able to distinguish the 
FVS and FVV curves from the FSS and FSV line of zero 
asymmetry. The SFF line peaks at below 10% making it 
also difficult to observe. 

The same plot for SPS 9 (figure [H bottom) shows low 
levels of asymmetry for all the curves, with the exception 
of the SFF curve. Its peak of over 10% asymmetry suggests 
that for these masses it could be picked out. It is unlikely 
that any of the other curves could be distinguished from 
each other for any of these mass spectra. 
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4 Model discrimination 

Here we apply the KuUback-Leibler distance 
notation, it is defined as 



In our 



KL(T, S) 



log 



p(to|T) 
p(to|S') 



p{m\T) dm, (5) 



where p{m\T) is the probabiUty density function for m 
given the distribution T, and analogously for p{m\S). The 
expression that distribution T is R times more likely than 
distribution S, on the basis of the data points {rui}, is 



P{T\{m..}) 
p{S\{mi})' 



(6) 



It was shown in [T^ that this can be rearranged to give 
a minimum number of events, N, needed such that dis- 
tribution T is calculated to be R times more likely than 
distribution S, assuming that T is the true distribution: 



N ^ 



logR 



KL(T, S) 



(7) 



in the limit A'' ^ 1 where p{S) and p{T) are the prior prob- 
abilities of each distribution. We must make an assump- 
tion about the true distribution as we must generate our 
data points for comparison from a particular distribution. 
This will not be the case when we have real data. We in- 
clude the factor \og{p{S)/p{T)) for completeness, however, 
we will set it to zero in our analysis. This is equivalent to 
assuming all distributions to be equally likely before we 
look at the data. Also, as pointed out in [1^, the result is 
invariant under diffeomorphisms m — *■ /(to), so the result 
would be unaffected if we had calculated with functions 
of for example, instead of to. 

The value N is an absolute lower bound on the required 
number of events. Once background and detector effects 
are included these will rise considerably, however, these 
effects vary from experiment to experiment and hence it 
is useful to have a universal lower bound. 

The results for the observable Pip distributions at SPS 
la (figure[5]) are given in tablelH The corresponding results 
for the curves at SPS 2 (figure El and SPS 9 (figure El) 
are shown in tables [5] and [SI The value R = 1000 has 
been taken here, so we are asking for one model to appear 
1000 times more likely than another. This corresponds to 
a 99.9% confidence level, but is a matter of choice. 

The lower numbers in the SPS 9 tables reflect the orig- 
inal impression from the graphs that the curves are easier 
to separate at this point than at SPS la or 2. The excep- 
tion is between the FSS and FSV curves, which was to 
be expected from the similar functional form. The values 
for SFF (which corresponds to the MSSM) are the lowest, 
but are still of the order of 100. These will be degraded in 
an experimental situation. 

The numbers in tables |4] ~ [6] treat the and W~ 
chains individually, however, we can reasonably expect 
that if one is observed, both will be. The relative numbers 
of the two chains again depend on the masses in the chain. 



With experimental data these values would be known, but 
here we rely on the MSSM values obtained in the same way 
as those in table[3] The results are shown in table[7l where 
the fraction of W chains which included a is denoted 
f±- 

When we consider both sets of data at once, equation 
(O is generalised to 



KL(T, S) 



-log 



p{m+\S+) 
p{m^\T^ 



p{m-\S- 



KL+(r, S) 



p{m \T ) Am 
KL-(T,5) (8) 



This gives the values of 



where p{m^\U^) = /± p{m\U^) 
N shown in table [51 

In order to illustrate how these numbers show an im- 
provement over treating the distributions separately, we 
consider the specific values of the (SFF,FVS) entry at 
SPS 9. For the P_ distribution alone it was 118, while 
for P+ alone it was 162. For this mass spectrum, a third 
of the chains have a . This means that looking at the 
whole sample, for every W~ event there are roughly two 
events. If the events contributed no discrimina- 
tory information (that is if ^(to+ISFF) =p(to+|FVS) for 



(a) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


697 


756 1237 


468 


FSS 


712 


oo 


177273 464 


1481 


FSV 


764 171080 


oo 498 


1687 


FVS 


1221 


434 


466 oo 


444 


fvv 


448 


1298 


1512 465 


oo 


(b) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


6728 


9459 975 


2801 


FSS 


7728 


oo 


177273 732 


1689 


FSV 


10408 171080 


oo 819 


2022 


FVS 


938 


688 


778 oo 


5523 


FVV 


2734 


1590 


1932 5605 


oo 



Table 4. Number of events needed, with SPS la masses, to 
disfavour the column model with respect to the row model by 
a factor of 1/1000, assuming the data to come from the row 
model, for (a) dP-/dm and (b) dP+/dm. 



(a) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


1220 


2223 704 


2166 


FSS 


1608 


oo 


19314 570 


1292 


FSV 


2668 17780 


oo 738 


2047 


FVS 


721 


560 


730 oo 


3181 


FVV 


2267 


1240 


2016 3211 


oo 


(b) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


1484 


1468 586 


649 


FSS 


1531 


oo 


19314 639 


1106 


FSV 


1483 17780 


oo 853 


1655 


FVS 


572 


619 


840 oo 


5551 


FVV 


630 


1081 


1638 5660 


oo 



Table 5. As in table |4l for the SPS 2 mass spectrum, for 
dP-/dm and (b) dP+/dm. 



8 



Jennifer M. Smillie: Spin Correlations in Decay Chains Involving W Bosons 



(a) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


90 


90 118 


87 


FSS 


83 


oo 5939353 648 


1686 


FSV 


83 5888890 


oo 659 


1734 


FVS 


97 


608 


618 oo 


1780 


FVV 


73 


1605 


1654 1844 


oo 


(b) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


123 


124 162 


121 


FSS 


117 


oo 5939353 666 


1686 


FSV 


118 5888890 


oo 677 


1735 


FVS 


139 


626 


637 oo 


2176 


FVV 


105 


1609 


1659 2253 


oo 



Table 6. As in table |4] for the SPS 9 mass spectrum, for (a) 
dP-/dm and (b) dP+/dm. 



Spectrum 


/+ 


f- 


SPS la 
SPS 2 
SPS 9 


0.57 
0.68 
0.67 


0.43 
0.32 
0.33 



Table 7. Fractions, f±, of total number of W chains which 
include a for each mass spectrum. 



(a) 


SFF 


FSS 


FSV FVS FVV 


SFF 


00 


1425 


1589 1073 


891 


FSS 


1476 


oo 


1.8-10^ 587 


1593 


FSV 


1619 1.8-10^ 


oo 642 


1863 


FVS 


1041 


549 


604 oo 


933 


FVV 


855 


1450 


1726 975 


oo 


(b) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


1388 


1647 619 


837 


FSS 


1554 


oo 


l.Q-W^ 615 


1160 


FSV 


1729 1.8-10" 


oo 812 


1763 


FVS 


613 


599 


801 oo 


4482 


FVV 


819 


1127 


1742 4550 


oo 


(c) 


SFF 


FSS 


FSV FVS FVV 


SFF 


oo 


110 


110 144 


107 


FSS 


103 


oo 5.9-10'^ 660 


1686 


FSV 


103 5.9-10'^ 


oo 671 


1734 


FVS 


122 


620 


631 oo 


2027 


FVV 


92 


1607 


1657 2100 


oo 



Table 8. Total number of and W~ events needed to dis- 
favour the column model with respect to the row model by 
factor of 1/1000, assuming data to come from the row model 
at (a) SPS la, (b) SPS 2 and (c) SPS 9. 



all TO"*"), then we would expect to need about 3 times the 
number of W~ events alone, 354. However, as the 
events do contribute to distinguishing the two models, we 
find that only 144 events in total are required. 

Equation ^ shows that when the prior probabilities of 
the models are equal (i.e. p{S) — p{T) as we have used) the 
number of events A^i calculated for a discrimination level 
i? = i?i is related to the number of events N2 calculated 
with i? = i?2 by a multiplicative factor. For example, to 
obtain the results for R = 20 (which corresponds to a 



95% confidence level), the numbers in table |4] should be 
multiplied by log(20)/ log(lOOO) ~ 0.43. 

It is instructive to consider how these events trans- 
late into required luminosity. The cross sections for these 
chains in the MSSM are given in tablelHl where the branch- 
ing ratios of X2 ~* ^ ^^'^ W ^ e,^ have been taken into 
account. This corresponds to considering the first row in 
each table — that in which the MSSM is the true sce- 
nario. The quoted required luminosity is calculated using 
the maximum number which appears in the first row of 
each table. 



Spectrum 


Cross Section (fb) 


Luminosity (fb ^) 


SPS la 


12.3 


129 


SPS 2 


1.41 


1171 


SPS 9 


0.03 


5473 



Table 9. Cross sections for chains of the form shown in figure 
[T]in the MSSM and corresponding integrated luminosity. 



The highest cross section is for the SPS la mass spec- 
trum, which is as expected as it is relatively light. This 
gives a required luminosity of 129 fb^^. The design in- 
tegrated luminosity for the LHC is 300 fb~^, before up- 
grade. This is encouraging, however, the required value 
will inevitably increase when detector and background ef- 
fects are considered. It looks unlikely that these studies 
could be conducted at this level of discrimination for ei- 
ther SPS 2 or SPS 9. The effect of the low numbers at 
SPS 9 has been suppressed by the small branching ratio 
of ql ^ <i'lX2 this point. 

5 Conclusions 

The spin correlations in the decay of a quark partner via 
a leptonic W boson decay, as exhibited in the invariant 
mass distributions of the quark and charged lepton, have 
been studied for three distinct SUSY-inspired mass spec- 
tra (SPS points la, 2 and 9). We have considered the 5 
possible spin assignments in the chain and studied the ex- 
tent to which they can be distinguished. The observable 
invariant mass distributions were constructed where we 
found that the distributions had similar functional form. 
The asymmetry constructed from these plots could be use- 
ful for distinguishing the SFF curve (which corresponds to 
the MSSM) from any of the other curves, but depends on 
the mass spectrum. 

The results were quantified using the KuUback-Leibler 
distance to give a lower bound on the number of events 
required to distinguish the spin assignments at a given 
level of certainty. This was applied to the distributions 
individually, and then to them combined. These provide 
a guide to how useful particular channels would be in a 
study like this. The lowest numbers were for the SPS 9 
mass spectrum where the lower bound was of the order of 
100 events when attempting to distinguish the SFF curve 
from others, and higher (in some cases considerably) when 
attempting to distinguish between the other distributions. 
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It therefore seems that this could be a useful method to 
distinguish the MSSM from other spin assignments in the 
chain, but will be less useful to distinguish amongst these 
alternatives. 

These bounds were converted to a luminosity require- 
ment for the case where the MSSM was the true scenario. 
The values for the SPS la mass spectrum were encourag- 
ing, while for SPS 2 and SPS 9 it appears unlikely that 
this method could give a high level of discrimination. 
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A Analytical formulae 

This section contains the formulae for the distributions 
plotted in section [3] (with the exception of FVV, see be- 
low). They are expressed in terms of two constants fci and 
^2 which are functions of the mass ratios described in sec- 
tion [31 



4y, 



ki = l + Y - Z and k2 ^ kl 
such that ki > k2 > 0. Then equations H]) and ^ 
ffii^ — fci (1 — cos 9 cos V') + k2 (cos 9 — cos ip) 



(9) 



2yY sin 9 sin cos < 



give 



(10) 



with maximum 2 (/ci -1-/02). We define the shorthands kf2 = 
ki ± k2 and fh = fhqi. 

Each distribution has different behaviour in the re- 
gions < fh? < 2k^2 ^-nd 2k^2 ^ ^ 2fcj*2, as can be 
seen in figures [2] - 14] and in the equations below. This is 
because high values of fh? can only occur for very specific 
angle configurations which cuts down the phase space and 
leads to logarithmic behaviour. This can be seen in the 
"No spin" curve in figures [2l^top),[3Utop) andjU^top) which 
represents the phase space. Without this effect, the "No 
Spin" distribution would continue linearly in m? > 2fc{^. 

SFF 

In the MSSM, the structure of the B-W-A vertex is 
l-|-a75 where a is defined by the parameters of the model. 
As this varies at each mass point, it is left explicitly in 
the equations below. Table [TUl lists the values of a at the 
particular points studied in this paper. 



Mass Spectrum 


SPS la 


SPS 2 


SPS 9 


a 


0.5083 


0.3875 


0.8155 



dP2 
dm 



Table 10. Numerical values of the parameter a for the differ- 
ent mass spectra studied in the text. 



16 {k2{{l + afki+A{l + afY 

+m2(l + Aa + a^ + (-1 + a^)^/Z)) 
-(2afciTO2 + (2(1 + a)2(l + fci) 

+ (l-|-a2)m2)y)iog(|±||)) 

< < 2k^2 

8fc?(l -I- af -I- (1 -I- (-6 -f a)a)rn^ 
x< +%ki{{l + af{k2+AY) (H) 
+n?{-2 - a{2 + a{2 + Vz)) + ^^Z)) 
+%ff?{-2(l + afY + k2{l + Aa + a^ 
-(1 - a^)Vz)) + 16y((l + af{b + 2fc2) 
-4(1 - a^)Vz) - m{2akim^ 
+ (2(l + a)2(l + fci) 
-f(l + a2)m2)y)iog(W^) 



2^12 ^ 



< 2k 



12 



3™ 



32A:2((1 + a2)(fc2 +2Y- SfciF) - 6Yy/Z{l - a^)) 
-16 (k2{m^ + 4y - fci(l - a)2 
+a(m2(4 + a)+ iaY) - n?y/Z{l - a^)) 
-{Y{2 + to2 - 4^/Z) + 2a{kifr? 2Y) 

+Y{2 + ff? + 4\/Z)a2) log 



{)<w? <2k 



12 



X < 



(12) 



8fc^(l - af + m4(l + a(6 + a)) 
-16Y{3 + 2k2- iVz 
+a(10 + a(3 + 2fc2 +4VZ))) 
+8ki{k2{l - af + 2m'^a - m'^VZ{l - a^) 
-4y(l + a^)) - 8m2(-2r(l + a^) 
+k2{l + 4a -f a2 + y/Z{-l + a^))) 
-|-16(y(2 + - 4\/Z) + 2a{kim^ + 2Y) 
+a^Y{2 + + iVZ)) log (^^^) 

2fcj"2 <fr? < 2k^2 



FSS 



dPi,2 
dm 



ifh 

2^2 



fcifc2 - 2yiog 



fcl+fc2 
ki — ^2 



< TO^ < 2fc 



12 



^{Qki - 2k2 - m2)(2(fci + ka) - ff?) ^^^^ 



16 V"'-! "'-^ 

-2riog(2(^ 



+ fc2) 
P 



2k^2 <ff?< '^k^2 



dPi 
dm 



3m 



32fc2((l + a^){kf + 2Y- 3kiY) - 6YVZ{1 - a^)) 
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FSV 

dPi,2 ^ 
dm 



3m 



2/c2(fc| + 12YZ) 

'fcifc2 + 2r(2Z-l)log(|±|j) 

< m^ < 2A:f2 

^ j^(6fci - 2A:2 - m2)(2(A:i + ^2) - m^) 
+2r(2Z-l)log(2i%^) 



(14) 



2fcf2 < m^ < 2k^^ 



FYS 

For the FVS chain the parameter a represents that in 
the C-q-B vertex, discussed at the end of section [51 




dA 

dm 



9m(fci + fe)^ 



8fc2(l + a2)(l + 2X){ki{ki + k2) - 2Y){kl + %Y) 
( 2fc2(l + a^)(fci +6m2) 

-4Xfc2m2(3 + a(3a- 2)) 
+ (m^(l + a2)(4fci + + X) 

-4r(l + a2-2X(l-a)2))log(f^) 

< m^ < 2A:f2 



Fig. 9. Distributions for the different scalar-vector-vector ver- 
tices described in the text. 



The new vertices give very similar distributions; the 
analytical expressions are given below. Due to the scalar 
in the chain the distributions for processes 1 and 2 are the 
same and therefore so are the distributions for the chains 
with positive and negative leptons. 



X < 



i ((2(2(fci + fcs) - m2)((l + a2)fc2(fci + k^) 

x(4fc2 +m^(15 - 16X)) 
-2(-(l + a2)(22fci + 26k2 + ISm^) 
+16((l-a)2(A:i +A:2) 
+ (l + a2)m2)X)r))/(fci + fc2)' 
+8((1 + a2)m2(4fci + m^){^l + X) 



(a) 



3m 



dm 



4fc2(fc| +6rz) 

' -fcifc2 + ((fci - 2y)2 + 2Y) log (fi±f|) 

< m^ < 2k^2 



X < 



-4r(i 



2X(l-a)2))log(H(^4+M)) 



2fci2 < < 2fcj*2 



^ (-(6fci - 2fc2 - m)(2(fci + fca) - m) 
+16((fci - 2y)2 + 2r) log (^^^^)) 



(17) 



dP2 
dm 



dm 



with 



(16) 



FVV 

The FVV distributions are too long to present here 
in a manageable way due to the complicated B^~W^- 
A vertex. They are available on request from the author. 
They also have the symmetry (ITS)) . 



dP, 



(b) 



dm 



3 m 



-fclfc2 + (fc2-2y)l0g(f^ 



+k2 



-k2 



< m^ < 2fc7 



(6fci 
x(2(fci 4 
+16{kl - 



- 2^2 — fh' 
^2) - m^) 
-2y)log( 

2k; 



2(fci+fc2) 



< m^ < 2fcV 



(18) 



B Higher dimensional couplings 

As mentioned in section [2l this analysis did not consider 
higher dimensional couplings with a different Lorentz struc- 
ture to that considered previously in the paper. For ex- 
ample, vertices with the form (a) p.qg^'^ — p'^q^ or (b) 
e^'^P'^Ppq^r are studied in the context of anomalous Higgs 
couplings [23] . 

Figure [9| shows the distributions for the FSV chain 
with these different vertices alongside the distribution shown 
before (marked FSV), for the SPS la mass spectrum. 



While the distributions are very similar to each other, 
they are quite different to that of the same chain with 
the lowest order vertex and to the distributions for the 
other chains shown in figure [5] Therefore vertices of this 
kind would be unlikely to be mistaken for those already 
discussed. 
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